Violation of time reversal symmetry in quantum noninvasive measurements 
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A quantum system in equilibrium is by definition in a state, which is invariant under time- 
reversal. However, in general projective quantum measurements of incompatible observables would 
obviously break the time-reversal symmetry, which one could try to circumvent by a weak quantum 
measurement leaving the state of the system intact. We show that, paradoxically, time-reversal 
symmetry is violated for this type of measurement, with the violation detectable in higher-order 
correlation functions. We propose an experiment with quantum dots to measure this apparent 
violation of time reversal symmetry in a third-order current correlation function. 



Microscopically, classical physics is symmetric under 
time reversalp] [2] — the arrow of time appears solely due 
to macroscopic loss of information which is the essence of 
the second law of thermodynamics. In quantum mechan- 
ics time reversal symmetry holds for standard measure- 
ments of compatible observables [3]. Strong, projective 
measurements of incompatible quantities affect future 
measurements and break the general time symmetry just 
as invasive classical measurements do, although a resid- 
ual symmetry (unrelated to the one discussed here) un- 
der exchange of pre- and postselected states survives [I]. 
However, the situation gets complicated when consider- 
ing noninvasive quantum measurements [6j. In this 
Letter we demonstrate that time reversal symmetry is 
broken in general for measurements of incompatible ob- 
servables even in the weak, non-disturbing limit. We also 
propose an experimentally feasible test of our conjecture 
in a quantum transport setup. Assuming experimen- 
tal confirmation, our result shows that quantum systems 
with time reversal symmetry will not appear to have time 
reversal symmetry. This is inconsistent with classical in- 
tuition and shows that the weak measurement limit in 
quantum mechanics cannot properly be considered non- 
invasive. 

Time reversal is often presented visually as reversing 
the direction of a moviepj. If the reversed film is com- 
patible with the equations of motion, the system is sym- 
metric under time reversal. The intuitiveness of this ex- 
planation is complicated by the necessity of making mea- 
surements to create the movie. Time reversal symmetry 
appears to be broken due to the second law of thermo- 
dynamics, which provides our arrow of time in the direc- 
tion of increasing entropy. This bias can be eliminated 
by considering stationary systems in thermal equilibrium 
where, for example, a backward-running clock would be 
just as likely as a forward-running clock. 

We consider an abstract "movie" to be the result 
of a series of measurements. Once the parameters of 
the system, measured quantities and times of measure- 
ments are all specified, the measurement results — the 
movies — are described by a joint probability density 
P(a\(tx), . . . , a n (t n )). A probabilistic description is re- 



quired by the statistical nature of the ensemble, mea- 
surement noise, and, in the case of quantum measure- 
ments, quantum randomness. The time reversed exper- 
iment involves reversed initial conditions (e.g. position 
r T = r and momentum p T = — p), external forces (e.g. 
magnetic field B T = — B), measured quantities a — > a T 
( T denotes the time reversal), and time — and therefore, 
ordering — of the measurements. Assuming for the mo- 
ment that the measurements have no effect on the time 
evolution, there will be an equivalence between the prob- 
ability densities of the forward (P) and reversed (P T ) 
sets of measurements: 

P(ai(ti),...,On(in)) =P T (<£(-i„),...,af(-ti)) . 

(1) 

This result leads to the principle of detailed balance [V, 
and reciprocity of thermodynamic fluxes [2J. Equation 
([I]) can be tested by comparing moments of the distri- 
bution. Of special interest is the moment where each 
measurement appears linearly: 

(ai(ti) • • • a n {t n )) = (c£(-tn) ■ ■ ■ al(-h)) T . (2) 

The time reversal operation is here applied to every mea- 
sured quantity, the parameters of the initial state and the 
evolution equations. In such a form, ([2| is independent of 
the validity of charge conjugation and parity symmetries 
and also of relativistic invariance [31 [H1[H]- For low-energy 
physics, time reversal alone is already a good symmetry 
in the equations of motion so, for equilibrium systems 
with no external magnetic field, P = P T . 

For equations ([!]) and ^ to hold, the measurements 
should not disturb the system, since any disturbance 
would create an asymmetry between before and after the 
measurement. On the other hand, quantum measure- 
ments of commuting observables clearly satisfy |T]), so 
non-disturbance of the system is too strong a require- 
ment: what matters is whether or not the measurements 
disturb each other. We adopt a straightforward defi- 
nition that a set of measurements is not disturbed by 
some additional operation or perturbation if the proba- 
bility density for the measurements is not affected. If a 
measurement can be added or removed without affecting 
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the joint probability density of the remaining measure- 
ments, then the complete set of measurements are mutu- 
ally non-disturbing. In other words, integrating over any 
single measurement should yield the same distribution 
that would be obtained if that measurement were never 
performed, 



J da,jP(ai, . . . ,a„) = P(ai 



tfj,...,a n ) (3) 



A set of measurements performed on an equilibrium 
ensemble and satisfying ([3| should therefore satisfy ([I]) 
and ([2]). This principle works not only for non- 
disturbing classical measurements, but also for quantum 
measurements of compatible observables (e.g. space- 
like separated[3])- These latter types of measurements 
disturb the system but are mutually non-disturbing in 
the sense of ([3]). Incompatible projective measurements 
[TU] do not satisfy ^ but also do not in general satisfy 
0. Nevertheless, for projective measurements one can 
find another symmetry — unrelated to ^ — under the ex- 
change of boundary conditions in pre- and postselected 
ensembles [U 

Weak measurements, introduced by Aharonov, Al- 
bert and Vaidman [5], however, allow for mutually non- 
disturbing measurements of incompatible quantities . 
They are also useful in refuting macrorealism by viola- 
tion of the Leggett-Garg inequality [IS], which has been 
recently tested experimentally [13] ■ In this scheme, the 
probability describing the measurement outcome P is a 
convolution P = D * Q of a large white noise D (in- 
dependent of the measured system, with uncertainty in- 
versely proportional to the strength of the measurement) 
with a quantum quasi-probability Q - a sometimes nega- 
tive Wigner-type distribution which is independent of the 
measurement strength |14j . The construction is based on 
positive operator- valued measures [15], using Gaussian 
Kraus operators [THl [T7] and taking the limit of infinite 
uncertainty. For the details, see [TS] and the Supple- 
mental Material A. Importantly, the pseudo-probability 
density for a series of weak measurements satisfies ([3]). 
Additionally, weak measurements result in only infinites- 
imal changes to the density matrix (and therefore, to the 
entropy) of any single system in the ensemble. Correla- 
tions of weak measurements are given by [151 EH] 

=Tr{A n (t n ),..{A 2 (t 2 ),{A 1 (t 1 ),p}}..}/2 n 

(4) 

for t n > ■ ■ ■ > t 2 > ti and {A, B} = AB + BA. Here 
A{t) denotes a quantum operator corresponding to the 
quantity a at time t and p is the initial density oper- 
ator. The striking paradox is that weak measurements 
satisfy |3]), at least as a limiting procedure, but it can 
easily be verified that weak measurement correlations Q 
in general violate the time reversal property ([2]). This 
asymmetry is only present for fast measurements of three 




FIG. 1: The double well, described effectively by two states 
|i) and jr), with energy shift 2e and tunneling amplitude r. 
Even in the ground state, quantum fluctuations allow jumps 
between the wells, which turn to be not time symmetric under 
noninvasive measurement. 



or more incompatible observables. In the case of compat- 
ible (AB = BA) or only two (not necessarily compatible) 
measurements the ordering is irrelevant and the symme- 
try (pi) holds. For slow measurements, each operator A(t) 
in BTis replaced with J f(t)A(t)dt, where f(t) turns on 
and off slowly compared to relevant timescales of the sys- 
tem. This slow measuring smoothes the resulting distri- 
bution Q so that any antisymmetric contributions vanish 
and therefore ([I]) and ^ will still apply. 

Let us demonstrate the paradox in a simple system 
consisting of a particle in a double-well potential as in 
Fig. [T] The particle is effectively described by the ground 
states of the left and right wells, \l) and |r) respectively. 
Higher excited states have much more energy and for 
low temperatures can be ignored, leaving an effective two 
state system. Using the basis states, the operator for 
expected location is Z — \l)(l\ — \r)(r\, and the effective 
Hamiltonian reads 



H = e(\l)(l\-\r)(r\) + T(\l)(r\ + \r}(l\), 



(5) 



where 2e is the energy difference between wells and r is 
the tunneling amplitude. Since no magnetic fields are 
present, H and Z must be even under time reversal. 
We are now in a position to test equation ([2| with z 
measured at three separate times and with the initial 
thermal state p oc exp(~H/kT). The correlation for 
three weak measurements can be calculated using Q and 
Z(t) = e im l h Ze~ lHt l h : 

(z(ti)z(t 2 )z(t 3 )) = a(e 2 + t 2 cos(2(i 3 - *a)A/ft)), (6) 



-(e/A 3 )tanh(A/A;T). For 



where A = V ' e 1 + r 2 , a 
this system and measurements, the expression corre- 
sponding to the right hand side of |2| differs from (6) 
only by the exchange of £3 — 1 2 with t 2 — 1\ . However, ( 6 ) 
is clearly asymmetric under this exchange, demonstrat- 
ing that time reversal symmetry is broken for correlations 
of quantum weak measurements. 

Despite the simplicity of the above example, a genuine, 
fast weak detection scheme is probably difficult to imple- 
ment experimentally in this case. Below, we present a 
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FIG. 2: (color online) (a) The system consisted of the dot 
(red) exchanging the elementary charge with a reservoir 
(blue), (b) The diagram of energy levels. At zero temper- 
ature, the Fermi sea of blue levels is full but electrons may 
still jump on and off the dot's red level. (c) Proposed detec- 
tion by an electric junction (yellow). The junction and dot 
are weakly coupled capacitively and the directly measured 
quantity is the current I through the junction biased with 
the voltage. The charge is allowed to jump between the dot 
and the lower reservoir but not the junction. 

more realistic example, leveraging recent developments in 
quantum dots |20) . We consider a quantum dot contain- 
ing a single energy level e, coupled to a Fermi reservoir 
by an energy independent coupling described effectively 
by the tunneling rate T/H, as depicted in Fig. [2^,,b. The 
occupation n on the dot (classically either or 1 in el- 
ementary charge units) is the measured observable N. 
The quantum observable and the Hamiltonian read [21] 

N = c ] c, H = eN + (7) 

J dE [^T/2nc i ^(E) + h.c. + E^(E)^{E)}, 

which describes energy-independent tunneling between 
the dot and reservoir, where e is dot level energy. We as- 
sume usual fermion anticommutation relations {tp, 4>} = 
0, {tf, 0} = if 1P ^ <j>, {c\c} = 1 and {^(E), iJj(E')} = 
S(E — E'). Spin is neglected here but if necessary all re- 
sults can be simply multiplied by 2. The initial state 
is p ex exp (-H/k B T). The Hamiltonian and the 
occupation are certainly symmetric under time reversal, 
H T = H and N T = N. 

To show the time asymmetry we will use the frequency 
domain, defining the third cumulant 

Sf {u), u)') = J dtdt' e iut+iu ' t ' ' (5n(t)Sn(t')6n(Q)) (8) 

with Sn — n — (n) . The asymmetry-probing quantity is 
the imaginary part of the third cumulant ImS^ (uj, to'), 
which should vanish if ^ holds. To calculate Q we use 



the close-time-path formalism [2"2l |2"5] , defining matrices 
in 2 x 2 Keldysh space 

*- (J !/«)•*-(£?) <f» 

with G R (cu) = ih/(hw - e + iT/2) = -G a *(lu) and 
G k {lu) = tanh(huj/2k B T)hT/(2(huj - e) 2 + T 2 /2). Then 

S»(uj,uj') = -J C ^TrG(a)Nx (10) 

[G(a + uj) + G(a + u}')}NG(a + uj + J)N. 

The integral can be performed analytically but the result 
contains digamma functions at finite temperatures. As 
suspected, ImS^ is not zero, see Fig. [3J Both imaginary 
and real parts vanish far from resonance. The asymme- 
try is the strongest at low temperatures (kT <§C e) and 
for comparable energy, tunneling and frequency scales 
(e ~ r ~ hut). This suggests that it is zero-point fluc- 
tuations of the charge jumping on and off the dot that 
are responsible for the asymmetry. The symmetry is re- 
stored if one of v, uj' , or w+w' is equal to 0. As expected, 
ImS 1 ^ vanishes for slow measurements us,uj' -^i T/H. In 
the limit uj,uj' —> 0, the result for S3 is a special case of 
the application of full counting statistics J23]. 

For the experimental confirmation of the asymmetry 
one must introduce a weakly coupled detector. We pro- 
pose an electric voltage-biased junction coupled weakly 
to the dot, so that its conductance depends on the charge 
on the dot, see Fig. |2j;. The externally measured quan- 
tity is the current / through the junction, in particular 
/ ~ 1° + where 1° is the intrinsic current in the junc- 
tion and x is its susceptibility due to the dot's charge. 
Then 5| ~ + X 3 Ss where is the internal current 
noise of the junction in the absence of the dot. Measure- 
ments of 5*1° have been demonstrated [24] therefore we 
expect the measurement of S3 to be feasible. Such an 
experiment will confirm the time reversal symmetry vi- 
olation only if the dot is not driven out of equilibrium. 
It is always possible for a certain parameter range, see 
Supplemental Material B, for the detailed model. 

We have shown that the property ([3| is not sufficient to 
guarantee time-reversal symmetries in sequences of quan- 
tum measurements. One could argue that the weak mea- 
surement still affects the system and forces time direction 
in this way. On the other hand, one expects a natural 
limit where the influence on the system is negligible and 
the time symmetry should hold. Although the entropy 
of the ensemble of many runs/copies of the systems is 
indeed changed by a finite value, even in the weak limit, 
we can safely assume the independence and identity be- 
tween systems so that the entropy change per system is 
small. 

In conclusion, we have shown that, contrary to classi- 
cal intuition, quantum noninvasive measurements violate 
time reversal symmetry. This violation is effectively a 
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n(u + w')/T e = r/2 
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FIG. 3: (color online) Time asymmetry of the third order cor- 



relation of dot occupation fluctuations S3 



in frequency 



domain at kT <C F. The symmetry is broken if Im S 3 7^ 0. 
(a) For e/F = 0.5 and arbitrary uj,ui' the asymmetry vanishes 
(white lines) at when oj, uj' or u + u' is equal to 0. The hexag- 
onal structure reflects the symmetry under permutations of 
frequencies, (b) Taking different values e/Y = 1,0.5,0.2 and 
uj — ui' it is clear that the maximal asymmetry occurs for 
energy /frequency parameters of the same order. 



failure of weak measurements to accurately reflect time- 
reversal symmetry inherent in a system. As such, it is 
independent of the validity of other symmetries such as 
charge-parity-time. Since quantum measurements of fi- 
nite strength manifestly break time reversal invariance, 
our result shows that, in contrast to classical measure- 
ments, all quantum measurements break time reversal 
invariance regardless of their strength. Finally, we em- 
phasize that only successful realization of the proposed 
or analogous experiment will ultimately confirm our con- 
jecture. 

We are grateful to Bertand Reulet, Christoph Bruder 
and Tomasz Dietl for helpful discussions. Financial sup- 
port from the DFG through SFB 767 and SP 1285 is 
acknowledged. 
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FIG. 4: The idea of weak measurement. The measured system 
(yellow) instantly interacts with a prepared ancilla (red wave- 
function), which is measured projectively afterwards (green 
detector). The procedure is repeated with identical but inde- 
pendent ancillas. 



A. WEAK MEASUREMENT 

To justify Eq. (4) in the manuscript we consider a se- 
ries of weak measurements. Following Aharonov, Albert 
and Vaidman [5], each weak measurement introduces an 
ancilla system and creates entanglement via an instant 
interaction Hamiltonian Hj = h6(t)gpA where g is the 
strength of interaction, p is momentum operator of the 
ancilla, conjugate to position x {[x,p\ = ih), and A is 
the measured observable. The interaction is followed by 
von Neumann projection [10] of the ancilla onto a posi- 
tion eigenstate which destroys the ancilla. The system 
can however be measured again with the next ancilla, as 
shown in Fig. [4J The density matrix after the 
surement is 



j th mea- 



Pj 



-mhh t 8 e iB & A i (A.l) 



where \<pj) is the initial prepared state of ancilla j. By 
inserting identity operations J da\a)(a\ — 1, the mea- 
surement interaction can be expressed as shifts of the 
ancilla wavefunction. 



Pj = / da 'j da"\<f>ji x j - 9ja'j))(<i>j(xj - 9ja")\ 

xi^x^i^xiom ( A - 2 ) 



In (A.2), the the state of ancilla j which has the shifted 
wavefunction 4>j(xj — gjd'j) is written as \<j)j(xj — gja'j)). 
The joint probability P(a\, . . . , a n ) =: P(a) is the proba- 
bility of measuring the ancillas in a set of position eigen- 
states with positions given by Xk = a-k/gk- This scaling 
makes the resulting distribution independent of measure- 



ment strength in both the strong and weak limits. 

P{a) = Tr|p„]^[^| afc / 5fc )(a fc / 5fc || (A.3) 
= J da' da" p n {a', a") 

(A.4) 

In ( |A.4[ ), pj is defined recursively by 
Pj{a[,a", d v a") = {a'^a'^-^pj-^a'-^a 1 -) . (A.5) 

Using Gaussian wavefunctions (f>{x) = (27r) -1 / 4 e~ x / 4 , a 
change of variables to a = (a' + a") /2 and da = a' — a" 
separates the joint probability density into a quasiprob- 
ability signal (Q) and detector noise (D). 



P{a)= J da D(a- a)Q(a), 



D(a - a) =Y[g k \0(gk( 



Q(a) = J dSa e -(9"5«) 2 /2 p n (g,,Sa) (A.6) 



The measure dda depends on a. Eq. (A.6 1 is the joint 



quasiprobability density for the series of von Neumann 
measurements. The quasiprobability has a well defined 
limit g — > 0. In this limit for time-resolved measure- 
ment, the averages with respect to this quasiprobability 
are given by 



(«!••• On) = 



I da'da"p(a',a")l[ a ^±^, (A.7) 



which is equivalent to (4). The genuine, measured prob- 
ability P — Q * D is positive definite because it contains 
also the large detection noise ~ l/g which is Gaussian, 
white and completely independent of the system. 

An alternative, equivalent approach is based on Gaus- 
sian POVMs and special Kraus operators [15-17]. Let us 
begin with the basic properties of POVM. The Kraus op- 
erators K(a) for an observable described by A with con- 
tinuous outcome a need only satisfy J daK^ (a)K(a) = 1. 
The act of measurement on the state defined by the den- 
sity matrix p results in the new state p(a) = K(a)pK'(a). 
The new state yields a normalized and positive definite 
probability density P(a) = Trp(a). The procedure can 
be repeated recursively for an arbitrary sequence of (not 
necessarily commuting) operators A±, . . . , A n , 

p{a\,...,a n ) = K(a n ) / 5(ai,...,a„_i)i^ t (a„) . (A. 8) 
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The corresponding probability density is given by 
P(ai, . . . ,a n ) = Tr p{a\, . . . , a n ). We now de- 
fine a family of Kraus operators, namely K\(a) — 
(2A/7r) 1 / 4 exp(-A(i - a) 2 ). It is clear that A — > oo 
should correspond to exact, strong, projective measure- 
ment, while A — > is a weak measurement and gives a 
large error. In fact, these Kraus operators are exactly 
those associated with the von Neumann measurements 
previously described (with A = g 2 ). We also see that 
strong projection changes the state (by collapse), while 
A — > gives /5(a) ~ p ; and hence this case corresponds 
to weak measurement. However, the repetition of the 
same measurement k times effectively means one mea- 
surement with A — > kX so, with k — > oo, even a weak 
coupling A -c 1 results in a strong measurement. For 
an arbitrary sequence of measurements, we can write the 
final density matrix as the convolution 



voltage 



Pa (a) = / da' g x (a')Y[gk(a k 

J l. 



(A.9) 



with g k {a) = e -2A fc a 2 ^/ 2 A fc /7r. Here A = (Ai,...,A„) 
a = (oi, . . . , a„), and da — da\ . . . da n . The quasi- 
density matrix g is given recursively by 



Qx(a) 



dt; 



,e-f / 2A " x 



2tT J y/2nX 



(A.10) 

i{t/2-j,)A n 



with the initial density matrix g — p for n — 0. We can 



interpret g in (A.9) as some internal noise of the detec- 
tors which, in the limit A — > 0, should not influence the 
system. We define the quasiprobability [18] Q\ = Tr g\ 
and abbreviate Q = Qq. In this limit (A. 10) reduces to 



g(a) 



2n 



, a n -i)e 



(A.11) 

Note that Qo...o,a = Q, so the last measurement does not 
need to be weak (it can be even a projection). The av- 
erages with respect to Q ar e easil y calculated by means 
of the generating function (A. 11), e.g. (o)q = Tr Ap, 
(ab) Q = Tr{i,S}p/2 ; (abc) Q = Tr C{B, {A, p}}/4 for 
a = (a, b, c). As a straightforward generalization to con- 
tinuous measurement, we obtain 



(oi(ti) • • • a n (t n )) Q = 
Trplii^),!---^.!^ 



i),A n (t n )}- 



(A.12) 



for time ordered observables, t\ < < ■ ■ 
above presented ways agree in the limit A 



< U 

J2 



Both 



0. 




FIG. 5: The model of detecting the dot's charge. An electric 
junction contains another dot with effective occupation JV', 
coupled capacitively to the measured dot. The fluctuations 
of the current / in the junction biased by the voltage V depend 
on the dot's occupation with the proportionality constant 
X- 



in a broad level regime. The complete Hamiltonian, con- 
sisting of the dot part (7), and the junction part, reads 



H + e'N' + H V + e 2 NN'/C + eVQ L + J dE x 
J2 [y/T/2n£iP A (E) + h.c. + E^ A (E)^ A (E)}, 

(B.l) 



A=L,R 

N L =4>[(E)ME), N' = £d 



where Nl is the total number of elementary charges e in 
the left reservoir, C is the capacitance between the dot 
and the QPC, V ,e' denote effective tunneling rate and 
level energy of the QPC and V is the bias voltage. 

We measure current fluctuations in the junction, I(t), 
with the current in Heisenberg picture defined as I(t) = 
—edNL{t)/dt. Such fluctuations have already been mea- 
sured experimentally at low and high frequencies [24]. 
Most of fluctuations are just generated by the shot noise 
in the junction. Now, we consider a finite, but still very 
large capacitance. We expect a contribution from the sys- 
tem dot's charge fluctuation to of the order C~ 3 . We 
assume separation of the system's and detector's charac- 
teristic frequency scales, namely 



(r,£,fc s T)<ey<(r',£'), 



(B.2) 



which also includes the broad level approximation for the 
detector's dot. There exists a special parameter range, 



B. DETECTOR 



An effective model of weakly detecting the dot's charge 
using an electric junction is shown in Fig. [5] The junc- 
tion is treated as another dot between two reservoirs but 



T'C > eV ^ 



e 

rc 



(B.3) 



where the coupling is strong enough to extract informa- 
tion about N(t) which is not blurred by feedback and 
cross-correlation terms (left inequality) , but weak enough 
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not to drive the system dot out of equilibrium (right 
inequality). In this limit the dominating contributions 
to the detector current's third cumulant are given by 
^o + x 3 ^ with 

S™ = T(l-T)(l-2T)e 4 V/h, 

X = -e 2 d(I)/Cds', (B.4) 



where (I) = Te 2 V/h and effective transmission T = 
r' 2 /(e' 2 + r' 2 ). Although the ~ x 3 term in S£ is much 
smaller than the first one, other terms, corresponding to 
cross correlations and back action, arc negligible com- 
pared to the last term. 



